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NOMENCLATURE 


A vane set area, m 2 

a 1 >32.33,34 "effective" lift or drag curve slope 

c d> c s, dra g and lift constants, N-sec/m 2 

c vane chord, m 

d G global vane set drag, N 

d section drag, N/m 

E (v*) energy spectrum function (m/sec) 2 m 

k»kl reduced frequencies 

H total height of one vane, m 

L chracteristic length scale for turbulence, m 

l G global vane set lift, N 

1 section lift, N/m 

N number of vanes in vane set 

R correlation function 

r n velocity ratio (V/V n ) across vane set 

max 

S spectral function 

velocity ratio (V/Vq ) along vane span 

max 

T(kj,k2) extended Sears function 

fi turbulence intensity 

t time, sec 

u axial velocity, m/sec 

v free stream velocity, m/sec 

V 0 time averaged velocity of flow parallel to tunnel centerline, m/sec 





w 

X 

y 

z 

z 

8 

8 * 

n 

a 

G ij 

9 1 

y 

V 

V* 

V 1 > v 2 * v 3 
5 

P 

o 

t 

*2 

*3 

*ij 

*1J 


velocity perpendicular to vane set stagger line, m/sec 

lateral velocity, m/sec 

vertical (vane span) coordinate, m 

streamwise coordinate, m 

total distance along stagger line, m 

cross-stream coordinate, m 

angle between free stream velocity vector and line perpendicular to 
stagger line, rad 

angle between free stream and tunnel centerline, rad 
cross-stream separation distance, m 
viscous drag coefficient 

P P 

three-dimensional turbulent velocity power spectrum, (m/sec) m J 
stagger angle, rad 

wave number of cross-stream velocity fluctuations, m - ^ 

wave number of spanwise velocity fluctuations, m - ^ 

wave number magnitude, m _1 

wave number components, m~^ 

spanwise separation distance, m 

air density, kg/nr* 

root-mean-squared (rms) value 

time delay, sec 

spanwise correlation function, m 
vane to vane correlation function 

p 

one-dimensional turbulent velocity power spectrum, (m/sec) m 
two-dimensional turbulent velocity power spectrum, (m/sec)^ 


IV 



characteristic "break" frequency, rad/sec 


fi b 

w frequency of axial velocity fluctuations, rad/sec 

Subscripts 

Dq global drag 

d section drag 

Lq global lift 

l section lift 

m mean value 

u axial velocity 

w lateral velocity 

1 conditions upstream of vane set 

2 conditions downstream of vane set 


v 



SUMMARY 


An analytical method has been developed to estimate the unsteady aerodynamic 
forces caused by flow field turbulence on a wind tunnel turning-vane cascade system 
(vane set). This method approximates dynamic lift and drag by linearly perturbing 
the appropriate steady-state force equations, assuming that the dynamic loads are 
due only to free-stream turbulence and that this turbulence is homogeneous, iso- 
tropic, and Gaussian. Correlation and unsteady aerodynamic effects are also incor- 
porated into the analytical model. Using these assumptions, equations relating 
dynamic lift and drag to flow turbulence, mean velocity, and vane set geometry are 
derived. From these equations, estimates for the power spectra and rms (root-mean- 
squared value, a) loading of both lift and drag can be determined. 


INTRODUCTION 


NASA Ames Research Center currently has a modification project under way to 
expand the capabilities and improve the aerodynamic characteristics of its 40- by 
80-Foot Wind Tunnel. Various aspects of this project have been reported in earlier 
papers, with Corsiglia et al. (ref. 1), being the most recent. One modification 
already completed is the installation of a new drive system, which increased the 
maximum attainable speed in the existing closed circuit tunnel from 100 to 150 m/sec 
(200 to 300 knots). Also, a nonreturn leg with a 24- by 37-m (80- by 120-ft) test 
section has been added. This new tunnel, which shares the drive system with the 
closed circuit facility, is designed for a maximum test section velocity of 50 m/sec 
(100 knots). A plan view of this entire facility, called the National Full-Scale 
Aerodynamic Complex (NFAC), is shown in figure 1. 

Located inside the NFAC are eight sets of turning-vane cascades (vane sets). 
These vane sets, built to turn air efficiently around corners and provide acceptable 
flow quality throughout the wind tunnel circuit, were originally designed to with- 
stand the aerodynamic loads occurring at the maximum tunnel velocity of 100 m/sec. 

To ensure the structural integrity of these vane sets at the new design's maximum 
velocity, it became necessary to develop a theoretical method for estimating the 
aerodynamic loads. Two types of aerodynamic loads were required to test individual 
vane strength and vane superstructure strength. Local loads were the spanwise loads 
over one vane section used to determine the structural strength requirements of the 
individual vanes. Global loads were the net aerodynamic forces on an entire vane 
set used in the structural analysis of the vane set superstructure. 
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It is important to realize that the NFAC is not a completely new facility. As 
a result, very conservative loads estimates could cause the structural analysis to 
incorrectly indicate a need for the reinforcement or renovation of the existing 
structures. Thus it was necessary to provide accurate estimates of loads to ensure 
structural integrity and to minimize cost. Given this accuracy requirement, analyt- 
ical procedures for estimating both steady and unsteady vane set forces were 
derived. The method used in estimating the dynamic loads is presented in this 
report. Steady-state load estimation procedures will be published separately. 


DYNAMIC FORCE EQUATIONS 


The development of the dynamic loads estimation procedure began with the 
assumption that dynamic forces could be estimated by linearly perturbing the appro- 
priate steady state force equations. For a turning-vane cascade system with uniform 
inflow and outflow and no steady vertical velocity component, the global vane set 
lift and drag are (ref. 2) 


L G = oV z A ^ tan 8 1 * tan s 2^ 


D g = \ pV 2 A(tan 2 $ 2 - tan 2 8^ + | pV^A 


where 


( 1 ) 

( 2 ) 


V = V, cos 8 
Z 1 1 


(3) 


the area of the vane set is A, and is the viscous-loss coefficient. The sign 
conventions used for lift and drag (Lq, Dq), for the pertinent velocities <V Z , V^), 
and for the inflow and outflow angles 82) are shown in figure 2.^ ^ 


By dividing equations (1) and (2) by the total span of the vanes and then 
incorporating equation (3), the sectional lift and drag equations can be derived 

a = pV 2 | cos 2 S^tan & - tan 8 2 ) (4) 

d = I pV ^ N ( cos2 8 i tan2 s 2 " sin2 8 i + ^ 

where Z is the total distance along the vane set stagger line and Z/N is the 
distance between vanes. 
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Note that the direction of drag is opposite to that shown in reference 2. 
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Dynamic force equations can then be derived by linearly perturbing equa- 
tions (4) and (5) about the variables and 8 -j, with 82 assumed constant, 
analysis yields the following equations 


This 


51 = b 1 V 1 SV 1 + b 2 V } 6B 1 
m m 


6 d = b V 1 5V 1 + b^V* 58 1 
m m 


where the subscript m designates a mean (steady state) value, and 


( 6 ) 

(7) 


Z 2 / 

b.j = 2 p jy cos 8 1 I tan 8 ^ - tan s 2 

h m\ m 


) 


( 8 ) 


b 2 = 


■•§(■ 


2 2 

cos 8 ^ - sin 8 ^ + 2 cos 8 ^ sin 8 ^ tan & 2 

m m m m 


) 


(9) 


■»§(< 


2 Q . 2 n . 2 „ > 

cos 8 ^ tan 8 2 - sin 8 ^ .+ n £ 


m 


m 


( 10 ) 


l(-‘ 


bi. = p 77 -cos 8 . sin 8 , tan 8 ^ - cos 8 , sin 8 , 

4 N \ 1 1 2 1 1 

mm mm 


) 


(ID 


From figure 2, it follows that the axial and lateral inflow velocities (u,w) can be 
defined as 


u = cos 8 * 
w = V ^ sin 8 * 


( 12 ) 

(13) 


where 8 * = 8 -j - 0-|. Perturbing these equations yields 


6 u = 6 V. cos 8 * - w 68 , 

1 m m 1 

6 w = 6 V. sin 8 * + u 68 , 

i m m 1 


(14) 

(15) 


Solving these two equations for 681 and 6 V^ and substituting them into equa- 
tions ( 6 ) and ( 7 ), the fluctuating lift and drag can now be expressed as 


61 = C. 6 w + C„ 6 u 
1 1 
w u 


(16) 


6 d = C . 6 w + C , 6 u 
d d 

w u 


(17) 
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where 


C l = 2 P V a l 

w 

(18) 

C l -- 5 » W 0 oa 2 

(19) 

U 


c d = 5 pV o ca 3 

w 

(20) 

c d = 2 pV 0 ca 4 

(21) 

u 


* (b 1 w m * b 2 “./'(I pV 0 c ) 

(22) 

= <b 1 u m ' b 2" B )/ (? pV 0°) 

(23) 

= <b 3 w m * b 4 u m )/^ pV 0 oj 

(24) 

= (b 3 u m - °v) 

(25) 

= V- cos 0. cos 0* /cos 6, 

0 1 ml 

m 

(26) 

= V A cos 6 1 sin 6* /cos 8- 
0 1 m l 

m 

(27) 

= V 1 cos 0 1 /cos 0 1 

(28) 


m m 


UNSTEADY AERODYNAMICS AND THREE-DIMENSIONAL EFFECTS 


Equations (16) and (17) are estimates of unsteady sectional lift and drag using 
a quasi-steady airfoil approach and assuming that lift and drag respond instan- 
taneously to flow field changes. It is known, however, that the lift on an airfoil 
in a fluctuating flow field does not change instantaneously, but is dependent on the 
frequency of that fluctuation (ref. 3). This frequency dependency can be described 
by a transfer function relating lift to angle of attack ( 6 B -j ) . It is also known 
that a transfer function exists between the lift and the pulsating flow velocity, 

6V i (ref. 4). For this analysis, we used a single transfer function for both 
effects. It was also assumed that this lift transfer function could be applied to 
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dynamic drag. This latter assumption is reasonable for momentum drag, since it is 
dependent solely on angle of attack. The applicability of this transfer function to 
viscous drag ( n ^ » however, is questionable and thus it must be assumed that is 

small compared to momentum drag. Using these assumptions, equations (16) and (17) 
become 


61 


= C SL 5w T 
w 


c i Su T 
u 


(29) 


6d = C , 6w T + C. 6u T (30) 

d d 

w u 

where Filotas' approximation (ref. 5) was used for the extended Sears function 
(T). This approximation assumes a sinusoidally varying velocity gust engaging a 
wing of infinite aspect ratio (no wall effects). It incorporates unsteady effects 
(Sears function) and three-dimensional (3-D) effects due to the gust hitting the 
vane obliquely. Any possible cascade effects on T were ignored. The transfer 
function is defined as 


|T(k 1 ,k 2 ) 



(3D 


where the reduced frequencies are 


k 1 = (jjc/2V q 

(32) 

k 2 = vc/2 

(33) 

2 2 2 
k = k^ + k 2 

(34) 


The frequency of the axial velocity fluctuation is w and the wave number of the 
lateral velocity fluctuation is v. The effects of k^ and k 2 on T can be seen 
in figure 3. 

Equations (29) and (30) give the dynamic lift and drag loads at one fre- 
quency. Assuming that velocity fluctuations (turbulence) can be described as a 
linear superposition of sinusoidal components, it follows that lift and drag can be 
expressed in terms of spectra incorporating all frequencies. Through Fourier analy- 
sis, it can be shown that the power spectra of lift and drag become 

S 4 U,v,u) = Cf |T| 2 S (uj,v,u) + C 2 |T| 2 S (»,v,m) + ?C C £ |T| 2 «(S (a >,v,u)) (35) 

w u u w 

S d (a»,V,y) = C 2 |T| 2 S w (w,v,y) + C 2 |T| 2 S u (w,v,y) + 2G d C d |T| 2 «(S uw (a,,v,y) ) (36) 

w a u a u a w 
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where S,, is the lateral, S,, is the axial, and S,„, is the cross-turbulent veloc- 
ity spectra. The third independent variable in these equations (y) is defined as 
the wave number for the vertical velocity fluctuation. The form of these velocity 
spectra will be discussed later. 


GLOBAL DYNAMIC LOADS ON A VANE SET 


As mentioned previously, one goal of this analysis was to estimate the global 
dynamic loads on a vane set. One way of presenting these global loads is with one- 
dimensional power spectra and their associated rms integrated loads. For global 
lift, this one-dimensional spectrum can be defined as 


S r (w) 


fcJ" 


R, (t)e 


-itox 


dx 


( 37 ) 


where R, 
l G 


is the correlation function for global lift with time delay x. This 


correlation function is defined as 


R l (x) = E[L G (t)L G (t + x)] 

G 

where E[ ] is the expected value and is global lift. 


( 38 ) 


It is possible to get a more useful form for the global correlation function by 
making use of some definitions. First, global lift can be defined as the summation 
of lift of all the individual vanes, such that 


L G 



( 39 ) 


where L n is the lift on one vane. By substituting this expression into equa- 
tion (38) the global and one vane correlation functions are related by 


R, (t) 

l g 


N N 


■ZE 


R l ( t , n ) 
mn 


( 40 ) 


where R r is the one vane lift correlation function at time delay x and lateral 
L mn 

vane separation n = z n - z m . 

A second definition states that the lift on one vane is a function of section 
lift such that 
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£ dx 


(41) 


L = 
n 


( 


where H is equal to the total vane height. If follows from this relationship that 

f H 2 f H 1 

r l ( x , n ) =1 I R a C x , n , C)dx 1 dx 2 (42) 

mn J Q J 0 mn 


where R 0 is the section lift correlation function at time delay t, lateral vane 


“mn 


separation n, and vertical separation £ = x^ - X 2 « Using equation (29), R ^ 


mn 


subsequently becomes a function of the lateral and axial velocity correlations (the 
cross correlation of the velocities, R uw , is zero for isotropic turbulence). 


H„ -H, 


R L (x,n) = 
mn 


n 

•'o J o 


C- C | T | n ( t , n , C )dx dx 


£ £ 
w w 
m n 


1 ~~2 


-H, 


n 

•'o J 0 


C C T| R ( t , n , 5 )dx . dx 


£ £ 
u u 
m n 


1 2 


(43) 


Mote that unlike the global lift correlation function 



which is dependent 


only on time separation (t), the turbulence correlation functions are dependent on 
separations in three dimensions (x,n,5). The lateral-turbulent velocity correlation 
function ( R w ) can be defined as the Fourier transform of its respective power spec- 
trum, such that 





S (w*,v,u)e 
w 


iw*r+iv£+iun 


dw* dv dy 


(44) 


The axial correlation function, R u , can be defined similarly. 

Incorporating nonuniform velocity profiles into the calculations can be done 
easily at this point by including variables in equations ( 18)— (21 ) . For example 


C £ = 2 pV 0 ca iV (x) (45 

w max 
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where V n is the peak axial velocity, r is the velocity ratio (V/V n ) at 
u max u max 

each vane n, and s(x) is the velocity ratio along the vane height. For simplic- 
ity, it is assumed that the r n profile is similar for all vane heights, and the 
s(x) profile is similar for all vanes. 


By combining equations (40), (43), (44), and (45), equation (37) becomes 
N N -H, 


i r v' v' f 2 r 1 

= 2 ~, Z E I J S(X 1 


)s(x_)r r 
2 m n 


m=1 n=1 0 w 0 


- t w - 

jjj d»* dv d„ dx 


1 dx 2 dx (46) 


However, it can be shown that 


+ CD -f-CO 

Z- S^(ui*,v,y)e w dcu* dT = J" <S(cu - ti>*)S^(u>*, v,y )du>* 


= S^w.v.y) 


(47) 


where 6( ) is the Dirac delta function. Thus the one-dimensional global lift 
spectrum becomes 


S. (w) = 
L G 


+ 03 

Jl 


$ 2 ( dy 


(48) 


where 


® 2 ( v) 


H~ -H, 


n 

o o 


iv(x 1 -x 2 ) 

s(x 1 )s(x 2 )e dx 1 dx 2 


1/2 


f 


s(x)e lvX dx 


(49) 
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* 3 (iO - 


MM , . 

v— \ r-\ y(z -z ) 

y y r r e ra n 
/ j l_j m n 

n =1 m =1 


1/2 


5>„ 


l\iz r 


n =1 


( 50 ) 


and S 5 is given by equation (35) 


$2 and <t> ^ can be thought of as the spanwise and vane-to-vane correlation 
functions, respectively. Their effect is to attenuate the lift and drag forces when 
vane dimensions become large compared to the turbulent velocity scale. Alternately, 
in the limit as the scale of the turbulent velocity eddies goes to infinity (v and 
y tend to zero), a quasi-steady state situation occurs and $2 and $ 3 S° to their 
maximum values. The effect of $2 anc * ^3 as a function of nondimensional wave 
number can be seen in figures 4 and 5. 


The approach developed in this section is also applicable to global dynamic 
drag. Using this approach, the one-dimensional global drag spectrum can be shown as 


+ co 

(w) = J*J* $2( v )®3^)S d (a),v,y)dv dy 


(51) 


where S d is given by equation ( 36 ). 

It should be noted that the power spectra of equations (48) and (51) have been 
left in the two-sided format (-■» < u < ®). One-sided spectra (0 < u> < ®), used for 
many applications, can easily be derived from these simply by doubling the spectral 
values for positive frequencies. 

The squares of the rms lift and drag are subsequently determined by integrating 
their respective power spectra over the entire frequency range, such that 


+ 00 

\ = J \ 


(oj)dw 


(52) 


r 


S n (a))dw 
G 


(53) 
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LOCAL LOADS 


Local dynamic loads can be determined by incorporating two simplifying assump- 
tions into the global analysis. First, because local loads are for one vane only, 
the cross-stream velocity profile is not necessary and the normalized vane-to-vane 
correlation function (^/N) becomes unity for all values of wave number p (see 
fig. 5). Secondly, integration done in the vertical direction need only be done for 
the small vane section of interest and not for a whole vane. 


TURBULENCE SPECTRA 


In order to estimate the sectional dynamic lift and drag power spectra of 
equations (35) and (36), it was necessary to find analytic forms for the lateral, 
axial, and cross-velocity power spectra. For simplicity, a turbulence model was 
chosen that assumed the turbulent velocities were homogeneous, isotropic and 
Gaussian. It was assumed, however, that small differences between the axial and the 
lateral rms velocity components could be accounted for by adjusting these constants 
in the model. Using these assumptions, the 3-D power spectra for turbulence can be 
expressed in terms of a basic energy spectrum function E(v*) by the following 
equation (refs. 6 and 7) 



E(v») 




V . V . ) 

1 J 


(54) 


where E(v*) is a scalar function that describes the turbulent energy density as a 
function of the wave number magnitude v*, are the wave number magnitude compo- 
nents (i = 1 to 3), is the Kronecker delta, and 

v* = |\>| = ( Vl + v 2 + V 3 ) (55) 

Two-dimensional and one-dimensional spectra and the integrated rms can be 
derived as follows 


VvV = 


). .( V ) 

1J 1 


2 

°i 


j: 

Jl 

j: 


0 lj d “3 


*ij dv 2 


*n dv i 


(56) 


(57) 


(58) 
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For the axial and vertical directions important to our analysis, directions 1 and 3 
respectively, equation (54) becomes 

©l^v) = — (v 2 + v 3 ) (59) 

Qo 3 (v) = " -j p (v^ + v 2 ) (60) 

4lIV * 


° 1 3 (v) = ^ (-v^,) (61) 
4TTV* J 

0^ is the axial, 0 33 is the vertical, and 0i 3 is the cross spectrum of the 
turbulence. A form of the energy spectrum function that applies to wind tunnel 
turbulence (ref. 8) is 


2 4 


2 8L 


L v 


E(v*) = al f 


/ 1 r 2 2.3 
( 1 + L v*) 


( 62 ) 


where is the rms velocity fluctuation and L is the characteristic length 

scale for axial turbulence. Plugging this function into equations ( 59 ) - ( 6 1 ) and 
letting a vary between the directions yields 


'11 


2 2L' 

°1 2 
Tt 


2 2 

v 2 + ^3 

/ . .2 2.3 

(1 + L v#) 


(63) 


0 33 = ° 


2 2LT 

3 , 2 


2 2 
V 1 » V 2 

(1 + L 2 v 2 ) 3 


(64) 


2L 5 ( - v 1 v 3 } 

13 " CT 1 a 3 2 , .2 23 

TT ( 1 + L Vg) 


(65) 


Performing the integrations necessary to determine the one-dimensional velocity 
spectra of equation (57) yields 


2 L 

’11 = °1 * 


1 


P P 

(1 + L vp 


( 66 ) 
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2 , (1 + 3tA> 2 ) 

*33 = °3 2 H ~ .2 2.2 (67) 

( 1 + L \>^ ) 

<t> ^ o = 0 (68) 

Since the one-dimensional cross spectra in equation (68) is identically equal 
to zero, it follows that the cross-spectral components of equations (35) and (36) 
will not contribute to the load spectra in equations (48) and (51), implying that 
the loading due to axial and vertical turbulence are statistically independent. 

For this analysis, the wave number components are equivalent to 


tions (35) and (36) 



V 1 

= 

u/V 


(69) 


v 2 

= 

V 


(70) 


V 3 

= 

P 


(71) 

now be 

put 

in 

the 

functional 

form required for equa- 

S u (u!,v, 

,u) = 

1 

V 

°11 

(u/V, v,p) 

(72) 

S w (id,v, 

,p) = 

1 

V 

0 33 

(u>/V, v,p) 

(73) 


If a characteristic "break" frequency is then defined such that = V/L, these 
equations reduce to 


2 2 V 

S u (ui,v,y) = o u 


(vV) - (wV)‘ 


2 ,„2 , „,2 , ...2 ' 2.3 

it (fl. + (vV) + (pV) + w) 


(74) 


„ , .2 fV (vV) 2 ♦ 0 ) 2 

S w (w,u,p) - a w 2 2 2 , „,2 2.3 

it (P, + (vv) + (pV) + 0 ) ) 


(75) 


These final two equations for the turbulent velocity spectra can now be substituted 
into equations (35) and (36). Note that these spectra can be completely described 
with just three experimentally measured quantities, 1) the break frequency, 2) the 
rms velocity, and 3) the mean flow velocity. 
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PHYSICAL SIGNIFICANCE OF BREAK FREQUENCY 


The one-dimensional axial turbulent velocity power spectrum is derived by 
integrating equation (74) 


S (01) 
u 



S^(u),v,y)dv dy 




IT 



From this latter equation it Follows that 


(76) 


S U ( u, = O b ) 1 

S (w = 0) = 2 

u 


(77) 


Thus, the character istic "break" Frequency n b is the Frequency on the one- 
dimensional axial turbulence spectrum where the power is equal to one-halF oF the 
value at w = 0. 


RESULTS 


Utilizing equations (3D, (35), (36), (49), (50), (74), and (75), equa- 
tions (48) and (5D can be numerically integrated to give estimates For dynamic liFt 
and drag. The parameters whose values must be prescribed in order to integrate 
these equations, together with actual values For a typical case (vane set 1 oF the 
NFAC), are listed in tables 1 and 2. Some oF these parameters come directly From 
Full-scale vane set geometry (c,H,Z,N, 8 ^ ) , some From extrapolations oF Full-scale 
data ( B I , p ) , some From small-scale model data ( >1 ^ ’ ve l° c ity proFiles), and 
others From more complete data gathered at other large wind tunnel Facilities 
(fi b ,TI w ,TI u ) . JustiFications For the use oF these inputs and complete loads results 
For the entire NFAC are being published separately. Predicted global liFt and drag 
spectra and rms loads For vane set 1 are shown in Figures 6 and 7. 

The global loads estimates From this analysis include Four Factors not Found in 
simpler methods: T, tfg, $ 3 » ancl the velocity proFiles. It is thereFore oF some 

interest to determine how much eFFect each oF these Factors has on the Final 
results. By neglecting all oF them, the loads equations revert to their simplest 
Form: quasi-steady aerodynamics, perFectly correlated turbulence, and uniForm 

velocity proFiles. The global liFt and drag spectra For the case oF vane set 1 are 
shown in Figures 8 and 9. The individual eFFects oF T, $> 2 » and $ 3 can then 6e 
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seen by comparing the rms loading for this simple baseline case with results 
obtained by independently incorporating each factor into the loads equations. The 
results are presented in table 3. It is apparent that each one of these factors has 
a significant attenuating effect on the dynamic loads estimates. If these effects 
are combined, the lift and drag spectra of figures 10 and 11 can be calculated. If 
the appropriate velocity profiles are then included, the spectra of figures 6 and 7 
are once again derived. It should be noted that by incorporating all of these 
effects, the rms lift estimate has been reduced by a factor of 11 when compared to 
the baseline case (table 3). A similar attenuation factor of 8.7 is observed for 
rms drag. 

The input parameters required to estimate local loads are nearly identical to 
those of the global loads case. Simplifications do occur, however, since velocity 
profiles are neglected and integration is done over one small vane section. Exam- 
ples of lift and drag spectra and rms loading for a 3.47-m section of vane set 1 are 
shown in figures 12 and 13. 


CONCLUDING REMARKS 


An analytical method has been developed to estimate the unsteady aerodynamic 
forces caused by turbulence on turning vane cascades. It is a refinement of the 
quasi-steady, perfectly correlated approach and includes correlation and unsteady 
aerodynamic effects. The use of this method results in significant attenuation of 
dynamic loads estimates relative to quasi-steady predictions, especially at higher 
frequencies. A typical example from one vane set in the NFAC shows attenuation 
factors of 11 and 8.7 for global lift and drag, respectively. This attenuation 
could be important since overestimating the vane set loads could cause the subse- 
quent structural analysis to incorrectly indicate a need for renovation of existing 
structures. 

Although this method is undoubtedly more precise than the simple, quasi-steady, 
perfectly correlated approach, there has been no experimental confirmation of its 
accuracy. It is therefore recommended that predictions derived from this analysis 
be supported by experimental measurements. 
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TABLE 1.- GLOBAL LOADS INPUT FOR VANE SET 1 


Parameter 

Symbol 

Input value 

Vane chord 

c 

1 .83 m 

Length of one vane 

H 

20.91 m 

Distance along stagger line 

Z 

46.86 m 

Number of vanes in set 

N 

50 

Viscous drag, coefficient 

n Sl 

0. 1 14 

Inflow angle between free stream and centerline 

6 1 

s 2 

0.733 rad (42°) 

Outflow angle between free stream and centerline 

-0.855 rad (-49°) 

Stagger angle 

A 

0.785 rad (45°) 

Maximum time-averaged velocity 

°max 

75.74 m/sec 

Air density 

p 

1.172 kg/m 3 

"Break" frequency 

fi b 

3.0 Hz 

Lateral turbulence intensity 

TI 

w 

0.156 

Axial turbulence intensity 

TI 
1 L n 

0.115 
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TABLE 2.- GLOBAL LOADS INPUT FOR 
VANE SET 1— VELOCITY PROFILES 


Vane 

no. 

Cross-stream 
velocity ratio 

r n 

Vane 

no. 

Cross-stream 
velocity ratio 

r n 

1 

0.579 

26 

0.653 

2 

.650 

27 

.697 

3 

.721 

28 

.756 

4 

.792 

29 

.798 

5 

.862 

30 

.840 

6 

.922 

31 

.883 

7 

.953 

32 

.925 

8 

.980 

33 

.968 

9 

.988 

34 

.991 

10 

.995 

35 

.997 

11 

.995 

36 

1.000 

12 

.995 

37 

.999 

13 

.995 

38 

.997 

14 

.995 

39 

.994 

15 

.995 

40 

.992 

16 

.995 

41 

.987 

17 

.993 

42 

.976 

18 

.986 

43 

.965 

19 

.971 

44 

.936 

20 

.947- 

45 

.902 

21 

.906 

46 

.863 

22 

.857 

47 

.793 

23 

.809 

48 

.712 

24 

.760 

49 

.623 

25 

.700 

50 

.502 

Normalized vertical 

Vertical velocity 

distance, x/H 

ratio, s(x) 


0.000 

0.000 


.013 


.378 


.125 


.675 


.225 


.885 


.316 


.971 


.325 


.994 


.471 


1.000 


.604 


1 .000 


.694 


.974 


.788 


.861 


.863 


.715 


.944 


.607 


.991 


.484 
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TABLE 3.- ATTENUATION DUE TO T, * 2 , $ 3 , AND VELOCITY PROFILES 


Case number 

Global rms lift estimate 

Global rms drag estimate 

Baseline global 
rms lift estimate 

Baseline global 
rms drag estimate 

1 . Baseline estimate 

1.0 

1.0 

(neglect T, « 2 , * 3 * 



and velocity profiles) 



2. With T only 

.515 

.484 

3. With $2 only 

.467 

.469 

4. With $2 only 

.354 

.422 

5. With T, « 2 , $3 

.115 

. 147 

6 . With T, ® 2 , $ 3 , 

.091 

.114 

and profiles 
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k 2 = 0 (SEARS FUNCTION) 


k-j = 0 (3-D EFFECTS) 


0 1 2 3 

k 

Figure 3.- Extended Sears function vs. reduced frequency. 



1 2 

uH/2tt 


Spanwise correlation vs. nondimens ional wave number for a 
uniform velocity profile (s(x) = 1). 


Figure 5.- Vane-to-vane correlation vs. nondimensional wave number for a 

uniform velocity profile (r^ = 1). 
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Figure 6.- Global lift spectrum for vane set 1, o = 165,000 N. 
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Figure 11.- Global drag spectrum neglecting velocity profiles, 

a = 158,400 N. 
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Figure 12.- Local lift spectrum for 3.47 m section of vane set 1, 

a = 2960 N. 
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